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Early  time  structuring) 
Plasma  instabilities. 


Linear  theory  shows  that  when  a  plasma  is  created  with  a  drift  relative  to  an  ambient  plasma  acroes  a  magnetic^ 
field,  a  robust  electromagnetic  instability  is  excited,  if  the  drift  speed  exceeds  the  local  Alfven  velocity  Vo  >  K»- 
Simulations  of  this  ion-ion  filamentation  instability  in  a  homogeneous  system  show  that  large  amplitude  waves  with 
SB) Bo  >  1  develop  nonlinearly,  giving  rise  to  large  density  clumping  due  to  the  non  gyrotropic  nature  of  the  ion 
velocity  distribution.  In  this  report  hybrid  simulations  are  performed  to  examine  the  consequences  of  the  ion;ion 
filamentation  instability  for  H.A.N.E.  applications.  Uniform  periodic  simulations  are  considered  as  well  as  simulations 
consisting  of  a  finite  sized  debris  plasma  interacting  with  a  background  plasma.  In  the  later  case  the  filamentation 
instability  can  be  examined  in  the  realistic  situation  where  a  high  Mach  number  collisionless  shock  is  generated. 
Of  particular  interest  is  whether  the  ion-ion  filamentation  instability  forms  long  lasting  structure  (on  a  time  scale 
greater  than  the  shock  formation  time)  in  the  debris  piston.  .The  simulations  indicate  that  large  amplitude  density 
perturbations  with  6np  n,*  (D  represents  the  debris  plasma  and  A  represents  the  air)  generated  by  the  ion-ion 
filamentation  instability  at  early  time  do  persist  in  the  piston,  especially  at  the  debris-air  interface.  The  filamentation 
instability  may  be  responsible  for  some  of  the  debris  and  background  plasma  structure  associated  with  a  H.A.N.E., 
_ although  separating  out  effects  due  to  other  instabilities,  such  as  the  Alfv_enjon_cyclQtron  instability,  is  most  difficult.  < 
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Abstract 


Linear  theory  shows  that  when  a  plasma  is  created  with  a  drift  relative  to  an  ambient 
plasma  across  a  magnetic  field,  a  robust  electromagnetic  instability  is  excited,  if  the  drift 
speed  exceeds  the  local  Alfven  velocity  Vo  >  Va •  Simulations  of  this  ion- ion  filamentation 
instability  in  a  homogeneous  system  show  that  large  amplitude  waves  with  SB /Bo  >  1 
develop  nonlinearly,  giving  rise  to  large  density  clumping  due  to  the  non  gyrotropic  nature 
of  the  ion  velocity  distribution.  In  this  report  hybrid  simulations  are  performed  to  examine 
the  consequences  of  the  ion-ion  filamentation  instability  for  H.A.N.E.  applications.  Uni¬ 
form  periodic  simulations  are  considered  as  well  as  simulations  consisting  of  a  finite  sized 
debris  plasma  interacting  with  a  background  plasma.  In  the  later  case  the  filamentation  in¬ 
stability  can  be  examined  in  the  realistic  situation  where  a  high  Mach  number  collisionless 
shock  is  generated.  Of  particular  interest  is  whether  the  ion-ion  filamentation  instability 
forms  long  lasting  structure  (on  a  time  scale  greater  than  the  shock  formation  time)  in 
the  debris  piston.  The  simulations  indicate  that  large  amplitude  density  perturbations 
with  Srio  ^  tia  (D  represents  the  debris  plasma  and  A  represents  the  air)  generated  by 
the  ion-ion  filamentation  instability  at  early  time  do  persist  in  the  piston,  especially  at 
the  debris-air  interface.  The  filamentation  instability  may  be  responsible  for  some  of  the 
debris  and  background  plasma  structure  associated  with  a  H.A.N.E.,  although  separating 
out  effects  due  to  other  instabilities,  such  as  the  Alfven  ion  cyclotron  instability,  is  most 
difficult. 
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I.  Introduction 


The  late  time  plasma  structure  resulting  from  a  H.A.N.E.  can  be  detrimental  to  com¬ 
munications  and  sensors  of  various  kinds.  In  this  regard  it  is  important  to  be  able  to 
characterize  the  scale  sizes  of  the  density  fluctuations  (‘structure”)  in  the  late  time  limit. 
However,  in  order  to  do  this  it  is  necessary  to  know  the  process  involved  in  the  structure 
formation.  The  question  arises  as  to  whether  the  late  time  structure  is  related  only  to  pro¬ 
cesses  that  occur  late  in  time  or  whether  the  early  time  behavior  of  the  initial  debris/air 
interaction  plays  a  part  in  determining  the  late  time  plasma  characteristics.  Recently,  work 
at  Austin  Research  Associates  [Thompson  et  al.,  1988]  has  focused  attention  on  the  ion-ion 
filamentation  instability  (hereafter  sometimes  referred  to  as  the  IIF  instability)  as  an  early 
time  candidate  which  may  be  important  for  determining  late  time  plasma  structure. 


The  purpose  of  this  report  is  to  examine  the  ion-ion  filamentation  instability  using 
hybrid  particle  simulations.  It  is  important  to  keep  in  mind  that  the  ion-ion  filamentation 
instability  is  closely  related  to  the  Alfven  ion  cyclotron  (AIC)  instability  [Davidson  and 
Ogden,  1975,  Tanaka  et  al,  19S3,  Tanaka  et  al.  1985;  Ambrosiano  and  Brecht,  1987;  Otani, 
1988;  and  references  therein],  an  instability  that  has  been  extensively  studied  using  both 
theory  and  simulations.  The  free  energy  for  both  instabilities  is  an  anisotropy  in  the  ion 
energy  with  £',j_/£'l||  >  1  where  _L  (  ||  )  is  the  direction  perpendicular  (parallel)  to  the 
ambient  magnetic  field.  The  linear  theory  for  these  instabilities  will  be  briefly  reviewed 
in  Section  2.  Previous  simulations  of  the  ion-ion  filamentation  instability  [Thompson  et 
al..  1988]  and  the  Alfven  ion  cyclotron  instability  [Ambosiano  and  Brecht,  1987]  in  a 
homogeneous  system  show  that  large  amplitude  waves  with  SB/B0  >  1  develop  for  both 
instabilities.  Large  density  clumping  only  occurs  for  the  filamentation  instability  due  to 
the  nongyrotropic  nature  of  the  ion  distribution  function.  Previous  simulations  of  the  AIC 
instability  in  a  shock  environment  showed  the  development  of  both  large  amplitude  waves 
w’ith  6B/B0  >  1  and  density  perturbations  with  6n/no  >  1  [Thomas  and  Brecht,  198G; 
Winske  and  Quest  1988].  No  simulations  to  look  specifically  for  the  ion-ion  filamentation 
instability  in  a  shock  environment  have  been  carried  out.  Later,  we  will  comment  on  w'hy 
it  was  not  observed  in  previous  simulations  where  it  might  have  been  excited. 


This  report  addresses  the  behavior  of  the  ion-ion  filamentation  instability  in  a  uniform 
system,  and  more  relevant  to  the  H.A.N.E.  problem,  its  evolution  in  a  situation  where 
a  high  Mach  number  shock  is  generated.  As  in  homogeneous  configurations,  the  shock 
simulations  show  that  large  amplitude  density  fluctuations  are  produced  in  both  plasma 
species.  This  structure  persists  even  when  the  shock  front  is  considerably  ahead  of  the 
debris  piston,  suggesting  that  the  ion-ion  filamentation  may  be  responsible  for  some  of 
the  ionized  structure  associated  with  a  H.A.N.E.  in  both  the  debris  and  the  background 
plasma  which  remains  behind  the  front  of  the  debris  piston.  However,  as  will  be  showm 
later,  distinguishing  structure  produced  by  the  IIF  from  that  associated  with  the  AIC  is 
very  difficult. 
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The  rest  of  the  report  includes  a  brief  review  of  linear  theory,  a  short  description  of  the 
simulation  model,  and  a  discussion  of  the  simulation  results  for  the  ion-ion  filamentation 
instability  in  both  homogeneous  systems  and  in  a  shock  environment.  The  final  section 
will  list  the  salient  results  and  outstanding  questions. 

II.  Linear  Theory 


Linear  theory  for  the  AIC  instability  is  well  described  in  Davidson  and  Ogden  [1975] 

and  Ambrosiano  and  Brecht  [1987],  while  properties  of  the  IIF  instability  are  worked  out 

in  Thompson  et  al.  [1988].  As  we  shall  show,  for  the  most  unstable  case,  namely  for  k 
— + 

parallel  to  J3,  the  two  modes  axe  essentially  equivalent. 


For  k  j|  B,  the  linear  dispersion  equation  for  transverse  electromagnetic  waves  is 

J-Sk* J  fvkvl±aFi/?n  ±  ™iFi  =  0  (1) 


u>  —  ±  f 1j 


where  Fj,  u>j  =  (4 ne^nj/mj)1/2,  and  flj  =  tjB0jmjC  axe  the  velocity  distribution  function, 
the  plasma  frequency,  and  the  cyclotron  frequency  for  the  j-th  species  of  mass  m^,  charge 
ej,  and  density  nj.  Assuming  Fj  can  be  factored  into  a  uj.  part  and  a  Maxwellian  Uy  part 
(with  thermal  speed  vj),  (1)  can  be  written  as 


u r 


!*2 


(2) 


where  the  argument  of  the  plasma  dispersion  function  (Z)  is  (u>  ±  Qj)/kvj  and  %  = 
\rnjv2}  and  T±j  =  /  ^mjv'^FxjCpv.  Thus,  whether  the  perpendicular  distribution  is  a 
Maxwellian  F±j  ~  exp(— a  ring  distribution,  F±j  ~  S(v±  -  Fr)  or  a  beam, 

F±j  ~  <5(vj.i)6(vj_2  —  Vj,),  does  not  matter  when  k  ||  B  (the  most  unstable  case):  the 
perpendicular  velocity  distribution  enters  only  through  the  second  velocity  moment,  “T±" . 
It  should  also  be  noted  that  the  basic  instability  mechanism  does  not  depend  on  whether 
the  driving  anisotropy  is  in  the  ions  or  the  electrons,  or  even  if  there  is  an  ambient 

magnetic  field  present.  (Hereafter,  we  shall  use  7j_/T||  to  mean  T,j./Ti||-)  Under  different 
conditions,  the  instability  has  been  called  various  names  (e.g.,  electron  Weibel,  electron 
whistler,  ion  Weibel,  ion  whistler,  electromagnetic  ion  cyclotron,  Alfven  ion  cyclotron, 
etc.).  It  has  been  extensively  studied  theoretically  and  is  a  favorite  mode  for  testing 
electromagnetic  simulations  codes  [e.g.,  Morse  and  Nielson,  1971;  Brackbill  and  Forslund, 
1982;  Otani,  1988]. 


However,  we  will  specialize  to  the  H.A.N.E.  case  and  consider  debris  ions  streaming 
relative  to  the  air  across  a  magnetic  field  B0,  as  depicted  in  Figure  1.  We  know  this 
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DEBRIS 


Figure  1.  Schematic  of  the  simulations  showing  the  ion-ion  filamentation  instability  (IIF) 
developing  in  the  debris  and  the  Alfven  ion  cyclotron  instability  (AIC)  at  the  snock 
front. 
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situation  cam  be  unstable  to  various  fast  growing,  electrostatic  instabilities  [Lampe  et  al., 
1975],  but  when  the  debris  velocity  exceeds  the  Alfven  speed,  they  axe  either  stable  or 
not  very  efficient  at  coupling  the  streauns  together.  Insteaid,  a  collisionless  shock  forms 
in  front  of  the  debris  [e.g.,  Thomas  and  Brecht,  1986].  Because  the  Mach  number  >  3, 
the  shock  reflects  ions  and  the  AIC  instability  is  excited  as  the  reflected  ions  form  a  ring 
distribution  at  the  shock  front,  giving  rise  to  T±/T\\  1  [Wu  et  al.,  1984].  However, 

Thompson  et  al.  [1988]  have  shown  that  for  the  debris  streaming  relative  to  the  air,  the 
ion-ion  filaxnentation  instability  cam  also  be  excited  with  a  very  large  growth  rate,  7  Q,. 

Because  it  does  not  give  rise  to  momentum  coupling  (the  wavevectors  axe  primarily  along 
B ,  not  the  streaming  direction),  it  has  not  been  given  serious  attention  in  the  past.  As 
we  have  previously  discussed,  these  two  instabilities  are  not  separate  entities,  but  as  their 
sources  of  free  energy  axe  different,  we  shall  treat  them  as  such  in  the  subsequent  discussion. 
In  fact,  in  the  Appendix  we  show  that  the  two  instabilities  are  practically  indistinguishable, 
as  will  be  apparent  in  the  simulations  to  follow. 


We  briefly  review  two  interesting  aspects  of  these  instabilities.  First,  we  consider 
their  linear  properties.  Assuming  for  simplicity  that  the  debris  and  air  ions  have  the 
same  density,  mass,  and  fij  =  8nn}T\\j/Bl  =  1,  we  consider  numerical  solutions  of  (2) 
for  different  values  of  Tx/T\\  =  4, 16,36.  (For  the  IIF  this  corresponds  to  the  ions  having 
velocities  ±VJ,  =  2, 4, 6^,4  ( v A  =Alfven  speed  using  the  total  ion  density).  Figure  2  shows 
growth  rates  normalized  to  fiA  versus  wavenumber  kc/uA.  For  a  given  value  of  T±/T\\,  the 
growth  rate  peaks  for  kc/ujA  ~  1,  and  falls  to  zero  as  k  — ♦  0  and  for  k  large.  Decreasing 
Tl/T||  decreases  the  maximum  growth  rate  and  the  range  of  wavenumbers  over  which  the 
instability  is  excited  (cf,  Figure  2  of  Davidson  and  Ogden,  [1975]).  Furthermore,  it  can  be 
shown  [Davidson  and  Ogden,  1975]  that  at  maximum  growth: 


2 _ (£il±.)i/ 2 

1 2  T|, > 


[cf  Eq.  (207)  of  Thompson  et  al.,  (1988)]  and  for  /?,  ~  1, 


kc 

Ui 


~  hl±y/ 2 
“  4  V 


(3) 


(4) 


Second,  we  briefly  recall  the  nonlinear  properties  of  the  instability.  Consistent  with 
both  quasilinear  theory  and  simulations,  as  the  instability  evolves,  the  growing  unstable 
waves  pitch  angle  scatter  the  ions.  This  reduces  the  anisotropy  and  shifts  the  dominant 
wavenumber  for  growth  to  smaller  values  (cf,  Figure  2).  Lemons  et  al.  [1979]  show  from  a 
general  thermodynamic  argument  for  the  analogous  electron  instability  that  the  maximum 
magnetic  field  fluctuations,  SB2 ,  are  given  in  terms  of  the  initial  parameters: 


SB 

B0 


(5) 


5 


0 


5.0 


kc/CDj 


Figure  2.  Numerical  solutions  of  (2)  for  the  ion-ion  filamentation  instability  for  V^/v^  =2,3, 
and  4;  equivalently  for  the  Aifven  ion  cyclotron  instability  with  Tx/T\\  =  4,16,25 
(0i  =  1)  showing  growth  rates  -y/Q,  versus  wavenumber,  Jfcc/w,. 


which  is  well  confirmed  in  the  simulations  of  Ambrosiano  and  Brecht  [1987]  (Figure  7) 
and  those  to  follow  in  this  report.  The  theory  also  predicts  that  a  residual  anisotropy  will 
remain  in  a  system  of  finite  length  L,  where  kf  =  2it/ L  and 


~  1 


+  (h£)’ 

Uli 


(6) 


III.  Simulations 


A.  Preliminaries 


The  simulation  model  uses  massless  fluid  electrons  and  particle  ions,  assumes  quasineu¬ 
trality  and  ignores  the  transverse  component  of  the  displacement  current.  Thus,  the  model 
retains  the  physics  necessary  to  describe  low  frequency  electromagnetic  phenomena  to¬ 
gether  with  ion  kinetic  effects  [Brecht  and  Thomas,  1988].  The  simulations  described  in 
this  paper  use  an  algorithm  for  solving  the  electromagnetic  field  equations  given  in  Quest 
[1988].  Similar  simulations  have  also  been  performed  using  the  predictor-corrector  method 
of  Harned  [1982].  The  results  obtained  are  substantially  independent  of  the  field  algorithm 
used. 


A  variety  of  simulations  have  been  performed,  including  1-D  and  2-D  simulations  of 
a  periodic  homogeneous  system,  and  2-D  simulations  appropriate  for  the  generation  of  a 
shock  wave  for  H.A.N.E.  initial  conditions.  Many  parameters  were  varied  including  the 
ratio  of  the  debris  plasma  to  the  background  plasma,  the  thermal  speeds  of  the  various 
components,  the  lengths  of  the  simulation  region,  and  the  use  cf  different  initial  conditions 
for  generation  of  the  instability  in  a  shock  environment.  We  show  only  very  few  of  the 
results,  just  enough  to  emphasize  the  most  important  points. 


For  the  homogeneous  runs  the  magnetic  field  is  along  the  z  direction,  the  simulation 
coordinate  for  the  1-D  runs.  The  2-D  runs  are  in  the  z-x  plane.  The  initial  debris  beam 
velocity  is  along  x,  with  the  background  ions  at  rest.  The  density  ratio  is  tid/tia  =  9.  The 
system  length  is  calculated  in  terms  of  c/u>,,  where  the  total  ion  density  is  used  to  calculate 
the  ion  plasma  frequency.  Both  species  have  the  same  mass  and  initially  are  cold,  with 
essentially  no  thermal  velocity.  Typically,  50-100  cells  are  used  in  the  1-D  runs,  with  10000 
simulation  particles;  a  50  x  50  grid  is  used  for  the  2-D  cases  (e.g.,  Case  6)  with  200,000 
particles.  As  shown  in  Table  1,  typical  cell  sizes  range  from  0.2  -0.5  c/u>i  with  time  steps 
of  Q,At  ~  0.01. 
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Filamentation 

Simulations 

case 

type 

Lz{nz) 

Lz(nz) 

/?x  ] 

r.x/r.j 

SB/Bo 

Sn/rio 

1 

1-DP 

wmmm 

ESI 

10®  21 

1.6 

4.4 

2 

1-DP 

19 

106  10 

2.6 

3.9 

3 

1-DP 

1 

9.0 

106  3.6 

3.3 

3.7 

4 

1-DP 

1 

9.0 

30  2 

1.5 

1.6 

5 

1-DP 

1 

90 

200  3.8 

8.7 

2.8 

6 

2-DP 

90 

200  2.5 

20 

2.8 

7 

2-NDP 

■ 

15 

5.2 

8 

2-NDP 

■SI 

15 

2.4 

Table  1.  List  of  some  of  the  simulations  indicating  the  simulation  parameters  and  results. 
Cases  1-5  are  one  dimensional  periodic  simulations,  case  6  is  a  two  dimensional  periodic 
simulation,  and  cases  7  and  8  are  two  dimensional  nonperiodic  simulations.  Cases  1-3 
are  initialized  with  a  finite  perturbations  in  mode  three  of  the  system,  making  the  runs 
essentially  single  mode  simulations.  The  other  simulations  were  initialized  with  random 
perturbations.  The  density  ratio  of  the  debris  to  the  background  plasma  was  9  for  the 
periodic  runs.  Case  7  has  a  nonuniform  debris  piston  with  a  maximum  ratio  of  9  to  1 
relative  to  the  background  plasma  (45  to  1  in  Case  8).  The  lengths  are  expressed  in 
terms  of  cj u.’p,  calculated  using  the  total  ion  density  except  for  cases  7  and  8,  which  are 
calculated  in  terms  of  the  background  ion  density.  The  numbers  in  parentheses  indicate 
the  number  of  simulation  cells  used.  The  ambient  magnetic  field  is  in  the  z  direction  for 
all  cases  except  for  cases  7  and  8,  where  it  is  in  the  x  direction.  The  initial  streaming 
direction  is  in  the  x  direction  for  all  simulations  except  for  cases  7  and  8,  where  it  is  in  the 
z  direction.  (3±_  is  the  initial  perpendicular  thermal  Beta  of  the  plasma.  The  first  entry  in 
the  temperature  ratio  column  gives  the  initial  temperature  anisotropy  whereas  the  second 
entry  gives  the  calculated  anisotropy  at  the  end  of  the  simulation.  SB /Bo  and  Sn/n0  are 
the  peak  fluctuations  in  the  magnetic  field  strength  and  total  plasma  density  as  observed 
in  the  simulation,  (norepresents  the  maximum  initial  plasma  density).  The  values  of  these 
fluctuations  at  the  end  of  the  simulations  is  somewhat  reduced. 
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For  the  inhomogeneous  runs  in  two  dimensions  a  debris  slug  of  parabolic  shape  in  the 
streaming  (z)  direction  with  thickness  Ld  and  uniform  in  the  direction  (x)  along  the  ambi¬ 
ent  magnetic  field  interacts  with  a  uniform  air  background.  The  maximum  debris  density 
relative  to  the  air  is  nine,  and  quantities  are  normalized  in  terms  of  the  air  parameters. 
The  initial  streaming  velocity  of  the  debris  relative  to  the  air  is  10u^,  where  the  Alfven 
speed  is  computed  using  the  air  density.  For  the  two  runs  to  be  shown  later,  300,000 
simulation  particles  were  employed  with  the  system  size  of  100  x  5 c/u>,  with  A z  =  0.5cu;, 
and  Ax  =  O.lc/w,  (to  emphasize  the  short  wavelengths  that  are  excited  along  x)  and 
ft.Af  =  0.14. 


B.  Periodic  Simulations 


In  this  section  we  describe  the  periodic  simulations,  which  w-ere  carried  out  in  order 
to  examine  the  linear  and  nonlinear  properties  of  the  ion-ion  filamentation  instability  in 
a  homogeneous  system.  We  have  varied  the  system  length,  the  ratio  of  the  debris  plasma 
to  the  background  plasma,  the  thermal  speeds  of  the  plasmas,  and  whether  the  system  is 
1-D  or  2-D.  Table  1  summarizes  the  simulation  parameters  and  some  results  of  the  various 
cases  run.  A  few  salient  features  of  the  results  are  shown  in  Figures  3-7. 


Figure  3  shows  time  histories  of  the  fluctuating  magnetic  field  energy  and  the  debris 
and  background  ion  kinetic  energy  for  cases  1-4.  All  four  cases  have  the  same  0,±  =  9.  In 
the  first  three  cases  T,±/T, ||  is  essentially  infinite.  In  each  case  the  simulation  is  initialized 
with  a  perturbation  in  mode  m=3,  which  corresponds  to  a  different  wavelength  as  the 
system  size  is  changed  in  each  case.  The  magnetic  field  fluctuation  energy  histories  in  all 
four  cases  show  intervals  where  the  growth  is  exponential.  The  measured  growth  rates  of 
the  first  two  cases  are  nearly  identical  (7/fl,  =  2.2  ±0.2  and  1.9  ±0.1)  consistent  with  (3), 
i.e.,  the  maximum  growth  rate  depends  only  on  Consistent  with  (4),  the  mode  with  the 
iargest  growth  rate  is  actually  shorter  than  those  imposed  in  these  cases.  Similarly,  one  sees 
from  Figure  3  the  reduction  in  the  growth  rate  of  cases  3  (7/fl,  —  1  -4  ±  0- 1 )  and  4  (7/n,  — 
0.12  ±  0.02)  as  the  unstable  wave  is  made  longer  (case  3)  and  as  the  initial  anisotropy 
(case  4)  is  reduced  is  evident,  again  pointing  out  the  limited  usefulness  of  (3)  and  (4).  A 
numerical  solution  of  (2)  yields  growth  rates  fo*  the  four  cases  (7/fti  =  1.7, 1.1, 0.6, 0.4), 
which  are  in  fair  agreement  with  the  simulations.  A  difficulty  with  comparing  linear  theory 
to  the  simulations  stems  from  the  fact  that  the  theory  assumes  the  calculation  is  done  in 
the  center  of  mass  frame,  whereas  the  simulations  are  done  in  the  rest  frame  of  the  air 
ions,  in  order  to  compare  the  results  with  the  inhomogeneous  simulations  to  be  described 
later.  Further,  it  is  evident  from  Table  1  (righthand  side  of  column  6)  that  when  only 
shorter  wavelengths  are  allowed  in  the  system  (by  reducing  the  system  length),  the  final 
ion  temperature  anisotropy  is  much  larger  than  if  longer  wavelengths  are  also  allowed.  This 
is  true  even  though  the  maximum  perturbations  in  the  field  and  density  are  only  weakly 
dependent  upon  the  system  lengths,  and  is  consistent  with  the  (6).  Also  seen  in  Table  1 
from  the  small  values  of  6B/B0  and  6n/n0  is  the  quenching  effect  that  smaller  initial  ion 
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Figure  3.  Energy  partition  as  a  function  of  time.  Exponential  growth  in  the  magnetic  field 
energy  occurs  for  cases  (1)  (2)  and  (3).  Case  (4)  does  not  have  exponential  growth 
because  the  growth  rate  is  smaller  than  the  gyrofrequency  and  the  mode  polarization 
is  changing  quickly  as  a  function  of  time.  Debris  (beam)  and  air  (background)  energies 
change  out  of  phase  due  to  gyromotion. 
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temperature  anisotropies  have  on  the  instability,  consistent  with  (5).  Also  shown  in  Figure 
3  are  the  time  histories  of  the  kinetic  energies  for  both  ion  species,  denoted  in  the  figure 
as  the  beam  (i.e.,  debris)  and  background  (air).  The  energies  of  both  ion  components 
oscillate  with  a  period  of  2 n/Cl,,  which  evidently  corresponds  to  gyromotion  of  the  debris 
about  the  magnetic  field  and  the  out  of  phase  response  of  the  background  plasma.  The 
fact  that  the  kinetic  energies  oscillate  again  points  out  that  the  instability  does  not  give 
rise  to  much  momentum  or  energy  coupling  of  the  two  species. 


Figures  4  and  5  show  profiles  in  z  of  the  transverse  components  of  the  magnetic  field 
(By  and  Bz)  as  well  as  the  magnitude  of  B,  and  the  debris,  air  (no  and  n^)  and  total  ion 
density  at  two  times  (f l,t  =  3.2  and  4.9)  for  Case  2.  This  simulation,  as  well  as  cases  1 
and  3,  is  essentially  a  single  mode  simulation  in  which  mode  3  was  initially  excited  with 
a  finite  amplitude,  while  the  other  modes  were  not.  Growth  of  the  instability  is  observed 
for  times  on  the  order  of  the  ion  cyclotron  period,  and  the  polarization  of  the  magnetic 
field  perturbation  changes  in  time  approximately  in  the  same  direction  as  the  instantaneous 
relative  streaming  direction.  In  this  case,  Bz  was  initially  perturbed,  while  in  the  left  hand 
panel  in  Figure  4  By  is  dominant  (Cl,t  =  3.2)  and  at  a  later  time  Bz  is  again  dominant 
(Cl,t  =  4.9)  [right  hand  panel],  Thompson  et  al.  [1988]  show  the  magnetic  field  will  be 
polarized  primarily  in  the  relative  streaming  direction  as  long  as  7  >•  fl,.  However,  since 
the  debris  ions  are  rotating  in  the  magnetic  field,  the  streaming  direction  and  hence  the 
sense  of  polarization  change  on  a  time  scale  comparable  to  Cl,. 


Nonlinearly,  large  amplitude  density  perturbations  develop.  Initially,  the  density  fluc¬ 
tuations  of  the  two  ion  species  are  out  of  phase,  with  each  species  occupying  areas  va¬ 
cated  by  the  other  species,  so  that  the  total  ion  density  is  roughly  constant  in  space 
(X)n,  =constant).  This  is  expected  for  parallel  propagating  electromagnetic  waves:  the 
total  density  fluctuations  are  second  order  in  the  field  perturbations.  Later  in  time  dur¬ 
ing  the  nonlinear  stage,  the  total  ion  density  also  develops  significant  fluctuations  so  that 
Sn/n  ~  6B/B.  As  is  evident  from  the  figure,  the  density  variations  become  so  large  that 
there  are  regions  where  the  density  of  each  species  is  essentially  zero.  Note  that  since  this  is 
an  electromagnetic  instability,  such  large  density  perturbations  do  not  lead  to  saturation 
of  the  instability.  Compared  to  cases  shown  later,  both  the  magnetic  field  and  density 
fluctuations  exhibit  a  rather  odd  behavior.  Instead  of  the  more  usual  cascade  to  longer 
wavelengths  (as  discussed  below),  shorter  waves  appear  later  in  time.  This  occurs  because 
the  simulation  was  initialized  artificially  with  a  perturbation  of  fixed  wavelength  (mode  3 
in  this  case).  Although  this  mode  grows  to  large  amplitude,  shorter  waves  that  grow  faster 
are  still  excited.  Before  they  are  quenched  by  the  reduction  of  the  anisotropy,  they  grow 
large  enough  to  contribute  noticeably  to  the  spectrum.  The  fact  that  these  modes  appear 
on  the  edges  of  the  debris  cavities  may  also  suggest  some  sort  of  parametric  decay  process 
is  occurring  in  this  case. 


11 


Figure  4.  The  magnetic  field  components  Bx  and  Bv  and  the  total  magnetic  field  at 
n,t  =  3.2  Heft  side  )  and Q,t  =  4.9  (right  side)  for  case  2.  Since  the  growth  occurs  over 
a  time  scale  Q.t  ~  7r,  the  dominant  magnetic  field  component  also  changes  in  time.  In 
this  and  all  subsequent  plots  the  magnetic  field  is  normalized  relative  to  the  ambient 
magnetic  field  strength. 


Figure  5.  The  debris  density,  the  air  density,  and  the  total  density  corresponchng  to 
Figure  4.  The  two  ion  species  tend  to  mutually  exclude  each  other  in  space.  In  this 
and  all  subsequent  plots  the  density  is  normalized  relative  to  the  total  (debris  plus 
background)  initial  density. 


Case  5  displays  the  more  common  mode  cascade  from  short  wavelength  to  long  wave¬ 
length  perturbations.  The  simulation  was  initialized  with  random  noise,  so  the  fastest 
growing  wavelength  is  very  small  due  to  the  large  ion  temperature  anisotropy,  shorter 
than  is  resolvable  by  the  simulation.  Therefore,  early  in  time  the  dominant  wavelength  is 
the  shortest  mode  resolvable  by  the  grid  with  a  wavelength  of  approximately  six  or  seven 
cells.  The  mode  cascade  is  seen  both  in  the  plots  of  By  and  no  in  Figure  6.  The  pres¬ 
ence  of  the  mode  cascade  suggests  that  the  short  wavelength  modes  which  grow  early  in 
time  are  irrelevant  for  the  long  time  appearance  of  the  nonlinear  state.  Thus,  the  physical 
mechanism  controlling  the  short  wavelength  cutoff  need  not  be  known  in  order  to  under¬ 
stand  the  dominant  long  time,  long  wavelength  state  of  the  instability.  The  cascade  to 
longer  wavelengths  has  been  seen  in  previous  simulations  of  the  filamentation  instability 
[Thompson  et  al.,  1988]  and  the  AIC  instability  [Ambrosiano  and  Brecht,  1987],  and  is  a 
well  known  phenomena  for  such  electromagnetic  instabilities  [Morse  and  Nielson,  1971]. 
The  evolution  of  the  instability  is  basically  quasilineax.  The  shorter  wavelength  modes 
grow  more  rapidly,  according  to  linear  theory.  This  removes  only  a  small  portion  of  the 
anisotropy,  leading  to  the  stabilization  of  these  modes  and  the  continued  growth  of  longer 
wavelength  modes.  The  longest  modes  in  the  system  remain  unstable  to  infinitesimal  ion 
temperature  anisotropies  [cf  Eq.  (6)].  Therefore,  the  long  time  behavior  of  the  instability 
is  dominated  by  wavelengths  much  longer  than  those  corresponding  to  the  fastest  growing 
mode. 


Finally,  Figure  7  displays  features  of  a  two  dimensional  simulation  corresponding  to 
case  5,  again  showing  By  and  no  at  three  times.  Note  that  in  this  case  and  those  to 
follow  that  the  scales  change  from  panel  to  panel.  There  is  a  mode  cascade,  easier  to  see 
in  the  density  than  in  By  because  of  the  change  of  the  polarization  of  the  wave  discussed 
earlier.  At  the  end  of  the  simulation,  the  density  perturbations  are  more  aligned  along  the 
magnetic  field.  The  maximum  perturbation  amplitudes  are  roughly  the  same  as  for  the 
one  dimensional  simulation  (Table  1),  although  the  density  depletions  are  not  as  evident 
in  the  figure.  While  structure  exists  across  x,  the  dominant  changes  (cascade)  occur  along 
z.  This  is  true  because  the  maximum  growth  rates  are  for  those  modes  with  fcj.  k\\-  In 
the  homogeneous  periodic  model  there  is  little  difference  induced  by  dimensional  effects, 
i.e.,  by  changing  the  size  of  Lz. 


C.  Inhomogeneous  Simulations 


In  order  to  make  the  simulations  more  relevant  to  the  H.A.N.E.  case  we  consider  the 
interaction  of  a  finite  debris  piston  with  a  uniform  magnetized  background  plasma,  similar 
to  that  of  Thomas  and  Brecht  [1986].  Two  types  of  initial  conditions  have  been  examined. 
In  one  case  the  magnetic  field  is  completely  excluded  from  the  debris  piston;  this  is  a 
good  approximation  for  the  early  stages  of  an  essentially  collisionless  H.A.N.E.  Another 
type  of  initial  condition  is  to  assume  that  the  ambient  magnetic  field  is  imbedded  in  the 
debris  piston.  This  is  not  realistic  for  the  early  stages  of  a  H.A.N.E.  but  it  may  be  quite 
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reasonable  for  the  leading  edge  of  the  debris  piston  due  to  anomalous  resistivity  effects. 
There  is  relatively  little  difference  in  the  results,  since  the  dominant  physics  occurs  at  the 
interface  between  the  debris  piston  and  the  background  plasma.  Thus,  we  will  only  show 
results  for  the  case  where  the  magnetic  field  is  imbedded  in  the  piston  region  initially. 
Figure  8  displays  3-D  perspective  plots  of  the  debris  ion  density  and  the  magnetic  field 
component  By  for  three  times  in  the  simulation  of  Case  7.  At  early  time  (Cltf  =  0.6)  the 
ion-ion  filament ation  instability  develops  throughout  the  body  of  the  piston;  consistent 
with  linear  theory,  waves  grow  on  a  time  scale  much  faster  than  the  ion  gyroperiod  and 
therefore  much  faster  them  the  time  for  shock  formation.  Only  small  ripples  can  be  seen 
in  the  debris  density,  but  it  is  evident  from  the  By  plot  that  short  wavelength  fluctuations 
occur  throughout  the  debris  region.  (The  waves  upstream  in  the  air  represent  only  noise.) 


Later  on  (f lit  =  1.5  ,  still  before  shock  formation),  the  waves  continue  to  grow  and  on  a 
time  scale  comparable  to  that  of  the  homogeneous  simulations,  strong  perturbations  in  the 
magnetic  field  and  the  density  arise,  particularly  at  the  debris-air  interface.  The  density 
and  magnetic  field  perturbations  are  somewhat  reduced  from  their  values  in  the  periodic 
homogeneous  simulations  in  most  of  the  interaction  region,  with  the  waves  at  the  leading 
interface  between  the  debris  and  the  background  plasma  larger  than  in  the  homogeneous 
simulations.  These  differences  are  not  a  two  dimensional  effect,  as  we  have  seen  that  the 
two  dimensional  periodic  simulations  reveal  similar  perturbation  amplitudes  to  the  one 
dimensional  periodic  simulations.  However,  there  are  other  possible  contributing  factors. 
First,  the  parameters  of  the  interaction  are  comparable  only  for  the  peak  in  the  debris  shell, 
and  not  everywhere  throughout  the  shell,  leading  to  smaller  perturbations  in  the  interior  of 
the  debris  shell.  Second,  the  magnetic  field  is  being  swept  up  by  the  debris  piston  in  these 
inhomogeneous  simulations,  which  could  compress  and  amplify  the  waves  at  the  leading 
interface  between  the  debris  and  the  background  and  reduce  the  perturbations  in  the  body 
of  the  debris  shell.  Third,  the  waves  at  the  debris-air  interface  may  be  enhanced  due  to 
the  strong  gradients.  Similar  enhancements  are  seen  at  the  ramp  of  high  Mach  number 
quasiperpendicular  shocks  [Thomas  and  Brecht,  1986;  Winske  and  Quest,  1988]  and  are 
thought  to  be  due  to  coupling  to  other  modes,  such  as  drift  mirror  and  compressional  drift 
waves.  However,  this  enhancement  may  be  an  artifact  of  2-D  [Thomas,  1989].  Again,  we 
note  that  on  this  short  time  scale  (f2j t  ~  2),  little  momentum  and  energy  is  exchanged 
between  the  two  species. 


In  order  to  tranfer  momentum  and  energy  from  the  debris  to  the  air,  a  shock  wave 
forms  out  in  front  of  the  debris  piston,  as  is  seen  in  the  bottom  panels  (Qit  =  4.3). 
At  the  front  of  this  collisionless  shock  wave  are  large  amplitude  magnetic  and  density 
perturbations  produced  from  the  AIC  instability,  as  discussed  previously  by  Thomas  and 
Brecht  [1986]  and  Winske  and  Quest  [1988].  Although  the  calculation  has  not  been  carried 
out  very  far,  it  appears  that  the  shock  is  unaffected  by  the  development  of  structure  in 
the  debris  piston.  The  largest  amplitude  perturbations  are  maintained  at  the  interface 
between  the  debris  piston  and  the  background  plasma,  which  is  now  somewhat  removed 
from  the  shock  front.  Note  that  the  maximum  value  of  the  debris  density  is  greater  than 
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Figure  8.  Three-dimensional  perspective  plots  of  Bt  and  no  at  =  0.6, 1.5, 4.3  for 
a  nonperiodic  simulation  with  parameters  corresponding  to  Figures  6-7  and  debris 
thickness  Lq  =  50 c/w*.  Early  in  time  the  magnetic  field  is  polarized  in  the  direction 
of  relative  streaming.  Due  to  gyromotion,  the  By  component  quickly  becomes  com¬ 
parable.  Late  in  time  the  largest  magnetic  field  perturbations  are  at  the  shock  front 
of  the  shocked  background  plasma,  while  the  density  perturbations  are  largest  at  the 
debris-air  interface. 
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5,  while  the  initial  debris  density  is  0.9  and  the  air  density  is  0.1.  Thus,  6no  is  greater 
than  5  times  the  initial  debris  density  and  more  than  50  times  the  air  density.  And  as  in 
the  homogeneous  simulations,  6n/n  ~  6B/B  1. 


At  this  time  the  mode  cascade  to  longer  wavelengths  (in  x)  can  be  seen  in  the  body 
of  the  debris  plasma,  again  in  general  agreement  with  the  one  dimensional  and  two  di¬ 
mensional  periodic  simulations.  A  mode  cascade  also  occurs  at  the  interface  between  the 
debris  and  the  background  plasma,  leading  to  the  longest  wavelength  possible  at  the  end 
of  the  simulation  in  the  debris  plasma.  Note,  however,  that  now  the  debris-air  interface 
is  no  longer  sharp;  debris  ions  have  diffused  upstream  of  the  interface  and  exhibit  strong 
variations  in  x.  Because  these  debris  ions  have  leaked  out  to  where  they  can  interact  with 
waves  produced  at  the  shock,  and  because  the  IIF  and  the  AIC  instabilities  have  similar 
linear  and  nonlinear  properties,  the  origin  of  the  debris  structure  in  this  region  is  unclear. 
One  also  sees  in  the  By  plot  that  waves  occur  at  the  shock  front  and  that  the  wavelength 
is  not  the  longest  possible  for  the  simulation  (i.e.,  no  mode  cascade).  This  is  because  the 
ion  temperature  anisotropy  at  the  shock  which  drives  the  AIC  instability  is  continually 
being  maintained  by  reflected  ions. 


Other  two  dimensioned  nonperiodic  simulations  have  been  done  varying  the  initial 
conditions,  the  system  lengths,  the  Alfven  Mach  number,  the  thermal  speeds,  and  the 
thickness  of  the  debris  piston.  The  general  features  of  the  simulations  are  largely  indepen¬ 
dent  of  these  parameters  as  long  as  Ma  3,  the  directed  velocity  of  the  debris  is  much 
larger  than  the  thermal  speeds  of  the  ions,  and  the  system  length  includes  relevant  wave 
lengths.  Figure  9  shows  the  results  of  mother  run  (Case  8),  with  the  same  parameters  as 
Case  7  (same  debris  mass  and  velocity)  but  with  a  much  narrower  piston  initially  (lOc/u;,). 
Again,  fluctuations  due  to  the  IIF  are  produced  at  early  time  before  the  shock  forms,  with 
amplitudes  comparable  to  the  previous  case.  (The  waves  in  the  debris  look  smaller  because 
the  maximum  debris  density  is  much  larger).  At  later  times  there  are  again  three  regions: 
(1)  waves  at  the  shock  front  due  to  the  IAC,  (2)  waves  at  the  debris-air  interface,  which 
is  not  well  separated  yet  from  the  shock,  and  (3)  waves  of  smaller  amplitude  behind  the 
debris  slug  in  air  that  has  already  interacted  with  the  debris  via  the  IIF  instability  and 
has  been  left  behind.  At  later  times,  remnant  fluctuations  from  the  shock  will  continue  to 
interact  with  the  shocked  air,  so  that  it  is  not  possible  to  easily  determine  the  origin  of 
the  late  time  structure. 


Finally  we  note  that  the  presence  of  the  largest  amplitudes  at  the  interface  explains 
why  the  initial  conditions  regarding  the  magnetic  field  penetration  into  the  debris  piston  is 
largely  irrelevant  for  the  most  significant  structure.  For  realistic  H.A.N.E.  situations,  the 
interface  would  be  expected  to  be  the  most  important  for  the  early  time  structuring  since 
the  magnetic  field  will  not  have  penetrated  through  the  complete  debris  piston  early  in 
time.  We  again  remark  that  it  is  largely  irrelevant  whether  the  debris  is  being  controlled 
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by  one  instability  or  the  other;  the  important  point  is  that  the  piston  is  structured  on  ion 
inertial  length  scales  and  remains  so  on  an  interestingly  long  time  scale. 

IV.  Discussion 

This  report  has  demonstrated  several  properties  of  the  ion-ion  filamentation  instability. 
Overall,  the  instability  is  very  similar  to  the  AIC  instability  caused  by  1,  which 

is  not  surprising,  as  the  ion-ion  filamentation  instability  also  is  driven  by  an  anisotropy,  but 
due  to  relative  ion  streaming  across  a  magnetic  field.  Both  instabilities  have  similar  growth 
characteristics,  and  both  saturate  at  about  the  same  levels  in  terms  of  field  perturbations 
for  a  given  free  energy.  Both  instabilities  evolve  via  a  mode  cascade,  which  involves  short 
wavelength  perturbations  growing  first  and  removing  some  of  the  anisotropy  in  the  ion 
energy,  followed  by  slower  growing,  long  wavelength  perturbations  that  dominate  the  long 
time,  nonlinear  behavior.  It  is  during  this  later  stage  that  the  remaining  ion  energy 
anisotropy  is  removed. 

A  major  difference  in  the  two  instabilities  is  that  the  ion-ion  filamentation  instability 
produces  large  amplitude  density  perturbations  in  a  homogeneous  system  whereas  the  AIC 
instability  does  not.  This  is  apparently  due  to  the  fact  that  the  AIC  distribution  is  assumed 
to  be  gyrotropic  whereas  the  filamentation  distribution  function  is  not.  None  the  less,  it  is 
useful  to  remember  the  strong  connection  between  the  two  instabilities.  One  manifestation 
of  this  similarity  is  that  in  a  nonuniform  shock  system  both  instabilities  produce  significant 
density  structure  of  comparable  amplitude  and  wavelength. 

An  important  question  relating  to  H.A.N.E.’s  is  whether  the  structure  produced  by 
the  filamentation  instability  is  long-lived.  Periodic  simulations  show  that  the  density  struc¬ 
tures  cascade  to  long  wavelengths  where  they  remain  for  a  long  time.  However,  uniform 
periodic  simulations  are  of  questionable  value  for  this  particular  effect,  because  in  an  iso¬ 
lated  system  there  is  no  place  for  the  energy  to  go.  For  this  reason  two  dimensional 
simulations  appropriate  for  shock  formation  were  also  performed.  Simulations,  initialized 
assuming  either  that  the  magnetic  field  is  already  diffused  through  the  debris  plasma,  or 
alternatively  that  there  was  no  magnetic  field  in  the  debris  plasma  initially,  demonstrated 
rapid  growth  of  the  ion-ion  filamentation  instability  on  a  time  short  compared  to  shock 
formation.  The  initial  density  perturbations  persisted  in  a  compressed  and  amplified  form 
at  the  interface  between  the  debris  plasma  and  the  background  plasma,  even  after  the 
shock  wave  itself  was  located  far  ahead  of  the  debris  air  interface.  However,  because  of 
the  similar  characteristics  of  the  two  modes,  the  simulations  were  not  able  to  unambigu¬ 
ously  determine  whether  the  structure  in  the  debris  piston  was  associated  more  with  the 
ion-ion  filamentation  instability  which  occurred  before  shock  formation,  or  was  a  manifes¬ 
tation  of  the  AIC  instability  taking  place  in  the  disturbed  background  plasma  at  the  shock 
front.  In  addition,  the  largest  fluctuations  occur  in  the  region  of  steep  gradients  at  the 
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debris-air  interface,  suggesting  strong  coupling  to  other  modes,  such  as  the  drift  mirror 
and  compressional  drift  modes  [Migliuolo,  1983]. 


Other  simulations  of  debris-air  interactions  have  been  carried  out  in  recent  years,  but 
did  not  show  the  development  of  the  IIF  instability.  Goodrich  et  al.  [1985]  did  hybrid 
simulations  of  the  initial  debris-air  coupling,  but  because  the  calculations  were  1-D,  wTith 
B  (and  thus  k  for  the  IIF)  perpendicular  to  the  simulation  direction,  the  IIF  could  not 
have  developed.  Winske  [1986]  carried  out  particle  simulations  of  debris-air  interactions;  in 
most  of  the  runs  the  ions  were  unmagnetized,  so  again  the  IIF  could  not  be  excited.  Even 
in  cases  with  magnetized  ions,  however,  the  system  size,  chosen  to  resolve  short  wavelength 
electrostatic  modes,  was  probably  too  small  for  proper  growth  of  the  IIF.  Finally,  Thomas 
and  Brecht  [1986]  and  Winske  and  Quest  [1986]  did  2-D  hybrid  simulations,  similar  to  those 
described  in  Sec.  III.  In  many  of  these  cases  the  debris  shell  was  very  narrow  (~  lOc/u,',), 
and  rather  large  spatial  mesh  (A  ~  OAc/u;,)  and  time  steps  (fl,At  0.1)  were  used  to 
study  the  long  time  behavior  of  the  shock  that  forms  in  the  air.  These  again  contributed 
to  inhibit  the  growth  of  the  IIF.  While  large  fluctuations  were  observed  at  the  debris-air 
interface  on  the  time  scale  of  the  shock  formation,  they  were  attributed  to  the  AIC,  rather 
than  the  IIF  instability.  In  the  calculations  of  Thomas  and  Brecht  [1986]  any  development 
of  structure  in  the  debris  could  not  be  followed  in  time,  as  the  calculations  were  done  in 
the  shock  frame,  causing  the  debris  ions  to  eventually  leave  the  simulation  region. 


While  the  present  simulations  nicely  illustrate  structure  development  in  the  debris, 
we  should  also  point  out  some  limitations  of  the  calculations.  First,  we  have  considered 
only  cases  where  the  magnetic  field  is  exactly  perpendicular  to  the  streaming  direction. 
As  a  consequence,  the  debris  and  air  ions  stay  relatively  isolated,  which  allows  the  best 
separation  of  effects  due  to  the  IIF  and  the  AIC  instabilities.  In  the  more  realistic  case 
of  oblique  orientation  of  the  magnetic  field  relative  to  the  streaming  direction,  significant 
diffusion  of  the  debris  into  the  air  occurs.  Thus  makes  the  distinction  between  the  IIF 
and  the  AIC  nearly  impossible.  Second,  the  simulations  here  were  run  for  only  about 
one  gyroperiod,  just  long  enough  to  observe  the  growth  of  the  IIF  instability  and  the 
formation  of  the  hading  shock  wave.  Clearly,  to  investigate  the  persistence  of  the  debris 
structure  the  simulations  will  have  to  be  run  for  much  longer  periods  of  time.  Third,  the 
two  dimensional  nature  of  the  simulations  imposes  unphysical  constraints.  The  conserved 
canonical  momentum  in  the  symmetry  plane  of  the  simulations  prevents  ions  from  moving 
more  than  about  an  ion  gyroradius  from  the  burst  field  line,  and  so  there  is  no  possibility 
for  the  macroscopic  shock  behavior  to  be  affected  by  the  debris  clumps.  In  a  real  three 
dimensional  system,  this  is  no  longer  the  case  and  the  initial  debris  perturbations  may 
have  some  effect  on  the  shock  front.  For  instance,  the  shock  may  produce  jets  of  debris 
material  in  certain  directions,  depending  upon  the  angular  mass  and  velocity  distribution 
of  the  debris  piston.  The  debris  clumps  may  also  seed  low  Mach  number  structuring 
instabilities  for  those  cases  where  the  interaction  changes  quickly  from  high  Ma  to  low 
Ma  (e.g.,  upward  expansion  of  the  burst). 
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In  conclusion,  the  simulations  presented  in  this  paper  suggest  that  the  ion-ion  fila- 
mentation  instability  may  generate  some  of  the  structured  ionization  in  both  the  debris 
and  the  background  plasma.  Perturbations  in  the  debris  density  much  larger  than  the 
ambient  density  (friQ  n^)  have  been  observed.  However,  it  should  be  kept  in  mind 
that  the  largest  waves  occur  at  the  debris- air  interface,  where  they  are  enhanced  by  the 
strong  gradients  and  perhaps  by  the  computational  constraint  of  2-D  and/or  the  initial 
conditions.  The  long  time  behavior  should  be  dominated  by  waves  of  comparable  ampli¬ 
tude  and  wavelengths  on  the  order  of  k\\c/upi  ~  1  Three  dimensional  calculations  will  be 
required  to  determine  if  the  ion- ion  filamentation  instability  has  any  globed  effect  upon 
shock  stability,  and  2-D  calculations  rim  to  much  longer  in  time  are  needed  to  realistically 
determine  whether  this  early  time  structure  impacts  the  late  time  state. 
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Appendix:  Comparison  of  Ion-Ion  Filamentation  and  Alfven 
Ion  Cyclotron  Instabilities  for  Debris  Piston/Shock 

We  consider  the  situation  depicted  in  Figure  1  in  more  detail.  The  debris  is  character¬ 
ized  by  density  np,  mass  mp,  velocity  Vp  _L  B„ ,  while  the  corresponding  air  quantities  are 
nA,  rnA,  and  X'a  =  0.  In  the  regime  where  the  debris  overruns  the  air,  the  mean  density, 
mass,  and  velocity  are  thus: 

fi  =  nA  +  nD  (.4.1) 

m  =  ( nAmA  +  nDmD)/n  (-4.2) 

V  =  ( nDmoVo  -f  nAmAVA)/hm  (.4.3) 

Recalling  that  for  both  the  AIC  and  IIF  (Section  2): 

i  -  <!%>•" 


Wi  4l  Til 


(-4.5) 


In  the  debris  region.  T\\  —  ^ nmt,  ,  where  v,  ^thermal  velocity  (for  simplicity  the  same  for 
both  species),  while 


rr*  -  -  r  *2  ^ 

1  =  -nm\  f  —  - n 


-nAmA(VA  -  V')2  +  -npmp(Vp  -  V)i 


Normalizing  to  air  quantities,  XlA  =  eB0/mAc,  u:A  =  {Ai:nAe2/mA  )'/2  and  defining  the 
Mach  number  Mq  =  Vd/va  (where  v2A  —  Bl/4imAmA)  and  0A  =  Strn ATu/ B2.  one  finds 
for  the  IIF: 


7 

1  (1  +  AT)(ArA/)1/2 
(1  +  NM)3/2  *  D 

(.4.7 

fcc 

3  (1  -f  N)(N M)x!2  Md 

( -4  S 

4  (I  +  .VM)3/2  p'J2 

where  A  —  np/n^,  M  —  mo/TnA ■  Taking  N  1  and  M  ~  2,  this  reduces  to 


7  ^  Mp 

ft*  ~  2V2 

kc  ^  3  Mo 

uA  ~  8  ^y2 


(.4.9) 

(.4.10) 
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In  the  shock  region  the  debris  is  irrelevant;  the  AIC  is  driven  by  reflected  ions,  with 
density  nr  =  an  A  and  velocity  Vr  ~  Vd,  so  that 

nT±_  ~  -nrmAVl  ~  -nAmAV[,  (>1.11) 

Proceeding  as  before,  we  find  for  the  AIC 

J~  (A.  12) 


kc_  _  3 a1/2  Mp 

^  ~~  4 


(A.13) 


Because  for  high  Mach  number  interactions  a  — ►  1/3  [Quest,  1986],  it  is  evident  from 
comparing  (A. 9-10)  with  (A. 12-13)  that  the  growth  rates  and  wavenumbers  of  the  two 
instabilities  are  indistinguishable. 
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